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The branching fraction of the w” decays into charmless final states is estimated in the S- and 
D-wave charmonia mixing scheme. With the information of the hadronic decays of J/wW and w’, it 
is found that a large branching fraction, up to 13% of the total Y” decays may go to charmless final 
states. The experimental search for these decays is also discussed. 
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I. INTRODUCTION 


As more and more data are collected at 4)(3770) (short- 
ened as w”’), 27.7 pb—! by BES-II, 55 pb~! by CLEOc, 
and a few fb—' is in progress by CLEOc, the study of the 
non-DD decays of y” gets renaissance recently. 


Due to its mass coincides with the DD mass threshold, 
w" is believed to decay predominantly to DD, and the 
non-DD decays including transitions between charmo- 
nium states and decays into light hadrons are expected 
to be small. The first exclusive non-DD decay mode of 
w" was searched for using the Mark-ITI data sample [1], 
where evidence of Y” > J/yrntr was observed; after- 
wards theoretical calculation of its decay width was pub- 
lished [2]. This decay mode was searched for recently 
by BES and a signal of 3.lø was observed [3], which 
is in marginally agreement with the upper limit set by 
CLEOc [4]. The other searches of the non-DD decays 
of w" were all reported in either doctoral theses [1, 5] 
based on the Mark-III data or conference talk [6] based 
on the BES-II data, and no statistically significant re- 
sults were given. This may indicate that the data sam- 
ples utilized are still not large enough to search for these 
channels of small branching fractions. Now with much 
larger w"” samples, measurements of the cross section 
a(ete~ => y" + DD) = o(DD) have recently been 
reported by BES-II [6] and CLEOc [7] collaborations us- 
ing either single-tag or double-tag method. Comparing 
with the measured a(e+te~ — y” — anything), the dif- 
ference can be up to 1.4 nb (about 18% of the total cross 
section of the y” production in et e~ collision [8]), which 
indicates the existence of substantial non-DD decays for 
w". However, there is limited theoretical work on this 
aspect. In Ref. [8], it is estimated that at most 600 keV 
(~ 2.5%) of the y” total width of (23.6 +2.7) MeV is due 
to the radiative transition, and perhaps as much as an- 
other 100 keV (~ 0.4%) is due to the hadronic transition 
to J/prr. All these together are far from accounting for 
a deficit of 18% of the total y” width. 


In this paper, we concentrate on the charmless decays 
of". By charmless decay, we exclude those decay modes 


*Electronic address: wangp@mail.ihep.ac.cn 


with either open or hidden charm. We shall estimate 
the partial width of the charmless decay of w’' in the 
2S-1D charmonia mixing scenario, under the assumption 
that the perturbative QCD (pQCD) “12% rule” holds be- 
tween pure 1S and 25 charmonium states. The mixing 
between 2S and 1D states of charmonium was originally 
proposed to explain the measured large Fee of y”. It is 
suggested [9] that the mass eigenstates (3686) (short- 
ened as ~’) and =” are the mixtures of the 2S- and 1D- 
wave of charmonia, namely 7(2?.S,) and y(13D1) states. 
This was used by Rosner [10] in explaining the “pr puz- 
zle” in w' and J/w decays. He suggested that the mixing 
of y(23 S1) and y(13 D1) states is in such a way which 
leads to almost complete cancellation of the decay am- 
plitude of Y’ > pr, and the missing pr decay mode of 7’ 
shows up instead as enhanced decay mode of 7". This 
idea was then applied to solve the enhanced decay of 
w!' + K} K? relative to J/y > K2K9, and a prediction 
of the branching fraction of y” + KK? was given [11]. 
In principle, if this scenario is correct, it can be gener- 
alized and relate the partial widths of each individual 
mode in J/w, w' and ~" decays. Therefore by virtue of 
the branching fractions of J/q and 7)’ decays, we are able 
to estimate the corresponding decay branching fraction 
of =". This will be tested by the upcoming large y” data 
sample from CLEOc. 

In the following parts of the paper, we begin our study 
with a general review of the 12% rule between 7’ and 
J/y decays, followed by an introduction of the 2S- and 
1D-wave charmonia mixing scheme. Then we conduct 
the calculation on the possible charmless decay width of 
w" with the available information from J/w and 7’ de- 
cays. We show that the scenario and currently available 
information on J/~ and ~' decays can accommodate a 
partial width of 3.0 MeV charmless decays of y”. Finally 
we turn to the experimental searches for these decays. 
Both inclusive and exclusive methods are examined. 


I. 12% RULE AND 2S-1D MIXING 


From the pQCD, it is expected that both J/ and y' 
decaying into light hadrons are dominated by the annihi- 
lation of cē into three gluons, with widths proportional to 
the square of the wave function at the origin |(0)|? [12]. 


This yields the pQCD 12% rule, that is 
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The violation of the above rule was first observed in 
px and K**K7 + c.c. modes by Mark-II [13]. Since 
then BES-I, BES-II and CLEOc have measured many 
two-body decay modes of 7)’ [14-22]. Among them, some 
obey the rule , like baryon-antibaryon (BB) modes, while 
others are either suppressed in y’ decays, like vector- 
pseudoscalar (VP) and vector-tensor (VT) modes, or en- 
hanced, like K&K?. There have been many theoretical 
efforts trying to solve the puzzle [10, 23]. Among them, 
the 2.5-1D charmonia mixing scenario [10] predicts with 
little uncertainty B(Y” — pr) which agrees with experi- 
mental data [1, 24]. 

In S- and D-wave charmonia mixing scheme, the mass 
eigenstates 7’ and w are the mixtures of the S- and 
D-wave charmonia, namely 


IY’) = |23.5,) cos — |13 Dı) sin , 
ly”) = |2351) sin 0 + |18 D1) cosé , 


where 0 is the mixing angle between pure 7)(2°.S,) and 
w(1°D1) states and is fitted from the leptonic widths of 
w" and w’ to be either (—27 + 2)° or (12 + 2)° [10]. The 
latter value of 8 is consistent with the coupled channel 
estimates [9, 25] as well as the ratio between y’ and 1)” 
partial widths to J/yn*r [2]. Hereafter, the calcu- 
lations and discussions in this paper are solely for the 
mixing angle 0 = 12° [26]. 

As in the discussion of Ref. [10], since both hadronic 
and leptonic decay rates are proportional to the square 
of the wave function at the origin, it is expected that if 
w' is a pure 7)(2°S)) state, then for any hadronic final 
states f, 


T(w’ > ete) 
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(2) 
The electronic partial width of J/w is expressed in po- 

tential model by [27] 
4a? e? 


a a |Ras(0)}?, 


with a the QED fine structure constant, ee = 2/3, M jyy 
the J/y mass and Rıs(0) the radial 135S, wave function 
at the origin. 

Since y’ is not a pure (2381) state, its electronic par- 
tial width is expressed as [10] 
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with My the ~' mass, me the c-quark mass, R2s5(0) 
the radial 23, wave function at the origin and RYp(0) 


x |cos6R25(0) — 


the second derivative of the radial 183D, wave function 
at the origin. In the calculations in this paper, we 
take Rəs(0) = 0.734 GeV?/? and 5 Ri, (0)/(2V2m2) = 
0.095 GeV3/? from Refs. [10, 28]. 

If Eq. (2) holds for a pure 2S state, Y” > f, Y! > f 
and J/y —> f partial widths are to be [11] 
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M3, 
1 = Cy : 2 
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where C’y is a common factor for the final state f, My 
the y” mass, and 7 = |nje’® is a complex parameter 
with ¢ being the relative phase between (f|13 D1) and 
(f 2751). 

From Eq. (3), it is obvious that with Tr(J/y > f) 
and ['(v’ + f) known, two of the three parameters, Cy 
and complex 7, can be fixed, thus can be used to predict 
TQ" — f) with only one unknown parameter, say, the 
phase of 7. Thus the S- and D-wave mixing scenario 
provides a mathematical scheme to calculate the partial 
width of ù” decay to any exclusive final state, with its 
measured partial widths in J/y and w’ decays. 

However, the current information concerning the w! 
decay is extremely limited, which prevents us from esti- 
mating Qn values for most exclusive decay modes. Ta- 
ble I lists some hadronic final states which are measured 
both in J/w and y' decays, together with the calculated 
Qn defined in Eq. (1). Summing up all the channels in 
Table I makes less than 2% of the 7)’ decay through ggg 
annihilation. 

From Table I, we notice that compared with the 12% 
rule, the yw’ decays to 


1. the pseudoscalar-pseudoscalar (PP) mode K} K? is 
enhanced; 


2. the VP and VT modes are suppressed; 
3. most of the BB modes are consistent with it. 


The summed branching fractions and Qp values for these 
three categories of decay modes are evaluated and also 
listed in Table I. In estimating the charmless decays of 
4", we shall discuss these three different cases separately. 

Since the experimental information on the exclusive 
decays of 7’ is rather limited, we turn to inclusive 
branching fractions of J/y and ~' hadronic decays as 
an alternative. The estimation is based on the assump- 
tion that the decays of J/aw and %' in the lowest or- 
der of QCD are classified into hadronic decays (ggg), 
electromagnetic decays (7*), radiative decays into light 
hadrons (ygg), and transition to lower mass charmo- 
nium states (ccX) [31, 32]. Thus, using the relation 
B(ggg) + B(ygg) + B(y*) + B(cEX) = 1, one can derive 


TABLE I: Branching fractions and Qa values for some J/~ and ~(2S) decay channels. 


Modes Channels Byjy, (107?) By2s)(10~*) Qn (%) Ref. 
0707 TT 0.147 + 0.023 0.8 40.5 5435 [29] 
K+K- 0.237 + 0.031 1.0 0.7 42 + 30 [29] 
KLK? 0.182 + 0.014 0.52 + 0.07 28.8 + 3.7 [16, 17] 
sum PP 0.57 + 0.07 2.32 + 1.27 41.1 +22.8 
1707 pr 12.7 4 0.9 0.29 + 0.07 0.23 + 0.6 [30] 
K+K” (892)7 + c.c. 5.0 + 0.4 0.15 £0.08 0.340.2 
KOK” (892)? + cc. 4.2 +0.4 1.10 + 0.20 2.6 + 0.5 
wn? 0.42 + 0.06 0.20 + 0.06 48+1.6 
1727 w f2 (1270) 4.3 £0.6 2.05 + 0.56 4.8 1.5 [15] 
paz 10.9 + 2.2 2.55 + 0.87 2.31.0 
K* (892)? K3 (1430)? + c.c. 6.7 + 2.6 1.86 + 0.54 2.8 +1.3 
of} (1525) 1.23 +0.21 0.44 +0.16 3.6 + 1.4 
sum VP& VT 45.45 + 7.37 8.64 + 2.54 1.90 + 0.64 
BB pp 2.12 + 0.10 2.07 + 0.31 9841.5 [29] 
AX 1.30 + 0.12 1.81 + 0.34 13.9 + 2.9 
DE’ 1.27 + 0.17 1.2 40.6 9.4 4.9 
©(1385)*¥(1385)7 1.03 + 0.13 1.1+0.4 10.7 + 4.1 
BE 1.8 + 0.4 1.88 + 0.62} 10.4 4.2 
AtA 1.10 + 0.29 1.28 + 0.35 11.6 + 4.4 
sum BB 8.62 + 1.21 9.34 + 2.62 10.8 + 3.4 


Note: f simple normalization by EB = (1/2)B2. 


B(ggg) + B(ygg) by subtracting B(y*) and B(céX) from 
unity. 

The calculated values of B(y*) and B(céX), together 
with the values used to calculate them are summarized 
in Table II. As regards to =’, two final states yn(2S) and 
h-(1'P,)+X with faint branching fractions are neglected 
in our calculation. By deducting the contributions B(y*) 
and B(céX), we find that B(J/y > ggg) + B(J/ > 
799) = (73.5 £0.6)% and B(Y' > ggg) + Bly’ > 99) = 
(19.1 + 2.5)%, then the ratio of them is 


_ BY > ggg + 99) 


= ZY _ 7 IJI ee = = (26.0 + 3.5)% . 
B(J/b > gg9 +99) l M 


Qg (4) 


The above estimation is consistent with the previous 
ones [32, 34]. The relation between the decay rates of 
ggg and ygg is readily calculated in pQCD to the first 
order as [35] 


according to the recent results from BES [16, 17]. The 
relation in Eq. (4) was discussed in the literature as the 
hadronic excess in w’ decay [32, 34]. It implicates that 
while some modes are suppressed in y’ decays, the dom- 
inant part of ~' through ggg decays is enhanced relative 
to the 12% rule prediction in the light of J/w decays. 


TABLE II: Experimental data on the branching fractions for 
J/w and Y’ decays through virtual photon and to lower mass 
charmonium states used in this analysis. Most of the data are 
taken from PDG [29], except for B(J/w, Y —> 7* — hadrons), 
which are calculated by the product R- B(J/y, Y! > pty"), 
with R = 2.28 + 0.04 [33]. In estimating the errors of the 
sums, the correlations between the channels are considered. 


Channel BUMD Bw) 


x” — hadrons (13.4£0.33)%  (1.66-£0.18)% 


ete (5.9340.10)% (7.55+0.31) x 107? 
[(J/y > 799) _16 a (1 = 9%) utp (5.8840.10)% (7.3 +0.8)x107 
T(J/p > 999) 5 as(me) Had qtr” (2.8 +0.7) x1073 
75X (25.2240.43)%  (3.4340.27)% 
Using as(m-) = 0.28, one can estimate the ratio to be Ye (1.30.4)% (2.80.6) x107 
0.062. A similar relation can be deduced for the w’ de- nR Jf (31.7£1.1)% 
cays. So we obtain B(J/y —> ggg) ~ (69.2 £0.6)% and mm I /p (18.841.2)% 
Biy’ => ggg) ~ (18.04 2.4)%, while the “26.0% ratio” in dee eae an 
Eq. (4) stands well for both ggg and ygg. Although Q, mdf (9.6£2.1) x10 
is considerably enhanced relative to Qp in Eq. (1), it co- NEG (8.60.7)% 
å A A A 70 0 YXel (8.4+0.8)% 
incides with the ratio for the Ky.A; decay mode between ifs (6.4+0.6)% 
W' and J/w, which is aX (.340.4)% (77 442.5)% 


Ill. THE CHARMLESS DECAYS OF y” 


We define the enhancement or suppression factor as [8] 


1 +,- 
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In the 2S-1D mixing scheme, for any final state f, its 

partial width in y” decay can be related to its partial 

width in J/y and w' decays by Eq.(3), with an unknown 

parameter ¢ which is the phase of 7. This unknown phase 

constrains the predicted [(w" — f) in a finite range. We 

calculate 


Rp =T" > f)/TU/¥ > f) (7) 
as a function of Q(f) and plot it in Fig. 1. In the figure, 
the solid contour corresponds to the solution with ¢ = 0; 
the dashed one corresponds to the solution with @ = 
180°; and the hatched area corresponds to the solutions 
with ¢ taking other non-zero values. 

To make it clear, we discuss the final states in three 


situations: Q(f) < 1, Q(f) > 1, and Q(f) = 1. 


A. Final states with Q(f) <1 


If Q(f) < 1, the decay Y’ > f is suppressed relative 
to J/w —> f. The extreme situation is Q(f) — 0, cor- 
responding to the absence of the mode f in w’ decays. 
This is the case which was assumed for the pm mode in 
the original work to solve the pr puzzle by the S- and D- 
wave mixing [10]. If Q(f) = 0, the solution of the second 
equality of Eq. (3) simply yields 7 = R2s(0) cos6/sin@ 
which cannot have a non-zero phase, and Rp = 9.2. 

Generally, the suppression factor could be different 
from zero. Even the pa and other strongly suppressed 
VP modes are found in w’ decays recently by BESII [19] 
and CLEOc [22] with Q(V P) ~ O(10-7). In this case, 
there are two real and positive solutions of 7 as shown in 
Fig. 1 corresponding to the maximum and minimum of 
their possible partial widths in %” decays. The solutions 
with 7 having a non-zero phase yield the values of Rp 
between the minimum and maximum limits. 

For VT final states, which are measured to have 
Q(VT) ~ 1/3, with Eq. (3), we get 2.0 < Rr < 21.6, 
as shown in Fig. 1, where the upper and lower limits cor- 
respond to two real and positive solutions of 7, the range 
is due to the values of 7 with non-zero phases. 


B. Final states with Q(f) >1 


If Q(f) > 1, the decay w' > f is enhanced relative to 
J/y) — f. The extreme situation is Q(f) —> oo, corre- 
sponding to the complete absence of the final state f in 
J/ decays. For 
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FIG. 1: Rr versus Q(f). The solid contour corresponds to 
$ = 0; the dashed contour corresponds to ¢ = 180°; and the 
hatched area corresponds to ¢ having other non-zero values. 
The inset displays the variation of Rr in the vicinity of Q(f) = 
1, where the dot-dashed line denotes Rr = 1. 


there are two real solutions of 7, one is positive and the 
other negative. From the first equality of Eq. (3), it is 
seen that the positive solution leads to the larger Rp. i.e. 
larger T (Y” — f) (the solid contour in Fig. 1) while the 
negative solution leads to the smaller r(Y” — f) (the 
dashed contour in Fig. 1). 

For the known enhanced mode in y’ decays like K2K?, 
Q(K2K®) = 2.26, we find 1.4 < Rr < 52.5, which cor- 
responds to the %” decay partial width from 0.024 to 
0.87 keV. 

It should be noted that for finite T(w' > f), Q(f) > 
oo means the diminishing of [(J/w — f), which gives 
rise to Rr — oo according to the definition of Rp in 
Eq.(7). Under such circumstance, it is more intuitive to 
calculate 


Rp =r" > f)/T@' > f). 


Its variation as a function of Q(f) is shown in Fig. 2. As 
Q(f) > œ, the solid and dashed contours converge into 
the same point (Rp — 21). In such case S-wave state 
does not decay to f, but D-wave does. Its partial width 
in yw’ and w" decays comes solely from the contribution 
of 7, or the D-wave matrix element. 


C. Final states with Q(f) =1 


These final states observe the 12% rule in J/w and 
y' decays. In this case, apparently one solution is 
identical to the pure electromagnetic decays, with 7 = 
5Ri 5 (0) /(2V2m2) and Rr = 0.048. As in the leptonic 
decays, their partial widths in w’’ decays are small rela- 
tive to the partial widths in J/ decays. However, there 
are also other solutions with overwhelming contribution 
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FIG. 2: Rr versus Q(f). The solid contour corresponds to 
o = 0; dashed contour corresponds to ¢ = 180°; and the 
hatched area corresponds to ġ having other non-zero values. 


from 7 which lead to very large partial widths in y” de- 
cays. This can be seen from Fig. 1 at the point with 
Q(f) = 1. We get the maximum of Rp of 34.0 which 
corresponds to another real and positive value of 7. If 7 
has a non-zero phase, then 0.048 < Rp < 34.0. 


D. Numerical results 


From Fig. 1, we see that except in the range 0.52 < 
Q(f) < 2.06 and a small range of the phase, Rr is always 
greater than 1. This range excludes virtually all known 
decay modes except BB which has Q(BB) ~ 1. Even 
inside this range, there are other solutions by which Rp 
(or RE) is at O(10). It means in this scenario, contrary 
to the naive guess, the charmless decay width of 7)” is 
greater than that of J/y or w’. More surprising is that 
Rr could be as large as a few tens for the final states with 
Q(f) > 1. In general, the final states which are enhanced 
in 7’ decays possibly have a large combined partial width 
in y” decays, especially if the phase of 7 is zero or very 
small. 

There are reasons to assume that the phase ¢ which is 
between the matrix elements (f|2.S) and (f|1D) should 
be small [11]. For the decay mode like pz, since there 
is almost complete cancellation between cos @R25(0) and 
nsin@ so that (pm|w’) = cos@R25(0) — nsind x 0, the 
phase of 7 must be small. If this is to be extrapolated to 
all final states, the physics solution will follow the solid 
contour of Fig. 1. Another argument comes from the 
universal phase between the strong and electromagnetic 
amplitudes of the charmonium decays. It has been known 
that in the two-body decays of J/w, the phase between 
the strong and electromagnetic amplitudes is universally 
around 90° [31, 34]. Recently, it has been found that this 


phase is also consistent with the experimental data of 7’ 
and w" decays [24, 36]. Since there is no extra phase 
between 2.5 and 1D matrix elements due to electromag- 
netic interaction, as in the calculations of the leptonic 
decay rates of w’ and w", the universal phase between 
the strong and electromagnetic interactions implies there 
is no extra phase between the two matrix elements due 
to the strong interaction too, i.e. ¢ œ% 0. This conclusion 
means, for the modes which are enhanced in 7’ decays, 
their partial widths in y” decay must be greater than 
those in J/ or w' decays by more than an order of mag- 
nitude. 

In Sect. II, we estimate that B(J/w —> ggg) ~ (69.24 
0.6)% while B(Y' —> ggg) ~ (18.0 + 2.4)%. Among 
the final states, we know that VP and VT modes have 
Q(f) <1. For them, 


S B(J/) > VP, VT) 
S Biy > VP, VT) 


4.6% , 
8.6 x 1074. 


w 
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Furthermore, there are final states with Q(f) ~ 1 (such 


as BB), for them, 


> BUI / > BB) 
> BW! > BB) 


x 0.9%, 

x 9.3 x 1074. 

After subtracting the final states which are known to have 
Q(f) < 1 and Q(f) ~ 1, the remaining 63.8% of J/y 
decay with a total width of [(J/w —> r. f.s.) % 58.1 keV, 
and 17.8% of w’ decays with a total width of T(w' > 
r. f.s.) & 50.1 keV which are gluonic either has Q(f) > 1 
or Q(f) unknown. Here r.f.s. stands for the remaining 
final states. On the average, these final states have 


Q(r.f.s.) ~ 2.19 . 


This is roughly comparable to Q(K2K?) = 2.26. 

If there is no extra phase between 2S and 1D matrix 
elements as argued above, then Rp takes the maximum 
possible value with Rp x 51.6 for Q(r. f.s.) = 2.19. With 
this value, we find r(Y” —> r.f.s.) = Rr x T(J/y > 
r.f.s.) & 3.0 MeV for the partial width of these remaining 
final states in w’’ decays, which is 13% of the total w” 
width. 

The above calculation takes the averaged Q(f) for the 
final states with Q(f) > 1 and Q(f) unknown, so it 
merely serves as a rough estimation. The exact value 
of the partial width should be the sum of the individual 
final states which in general have various Q(f) values. 
At present a major impediment to do the accurate evalu- 
ation is the lack of experimental information. Neverthe- 
less, if we take 13% as charmless decays (the calculations 
done channel by channel for VP, VT and BB modes in 
Table I give a summed maximum possible width of 93 
keV in %” decays, or 0.4% of the y” total width) together 
with the charmonium transition contributions of 3% [8] 
(2.5% for radiative transition and 0.4% for maJ/w), we 
obtain a maximum of 7)" non-DD decay branching frac- 
tion of 16% in the 2S-1D mixing scenario to be compared 
with 18% as summarized in Ref. [8]. 


IV. EXPERIMENTAL MEASUREMENTS 


With the expected data of a few fb—! at y” from the 
running CLEOc and more from the future BES-III, it is 
of great importance to search for the charmless decays 
of w" experimentally. These measurements test the 2S- 
1D mixing scenario on the explanation of the pr puzzle 
in J/w and w' decays, and provide information on the 
relative phase between the 2S and 1D matrix elements, 
which is related to the charmonium decay dynamics. 

Besides the theoretical interests, such measurements 
are important for the experiment itself. So far, in the 
fitting of the y” resonance parameters with the scanned 
cross sections, it has been assumed that 7)’ decays com- 
pletely into DD [37-40]. If the non-DD branching frac- 
tion is substantial, it must be considered in the fitting in 
order to get self-consistent results. 


A. Inclusive method 


A substantial non-DD decays of Y”, including charm- 
less decays and charmonium transitions, originally 
caught attention from the comparison between the cross 
sections of the inclusive hadrons and DD at the y” peak. 
However, this is not unambiguous due to the poor statis- 
tics of the data samples and the complexity of the anal- 
ysis. In addition, results from different experiments are 
consistent with each other only marginally. 
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FIG. 3: The cross section o(e*e~ — hadrons) in the vicinity 


of the Y” resonance calculated with parameters provided by 
PDG. The total cross section a?! is divided into four parts: 
the non-resonance ow p, the radiative tails of J/w (oj;y) and 
Y (oy), and the Y” resonance (oyr). 


Fig. 3 shows o(ete~ — hadrons) in the vicinity of the 


w" resonance calculated with parameters by PDG [29]. 
The total cross section gt% can be expressed as 

or =ONRTOJ/y + Oy + Oy", (8) 
which contains four parts: the non-resonance cross sec- 
tion owr, the radiative tails of J/y (ozjy) and y' (oy), 
and the w” resonance (gy). The non-resonance cross 
section is usually expressed in terms of R. value and the 
u pair cross section as oyr = R-a(ete” — pty). The 
Breit-Wigner formula is adopted to depict resonances of 
J/, w' and w", where the total decay width of 7” is 
energy dependent: 
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where p is the D? or D+ momentum, r is the classical 
interaction radius, and Cr is defined as follows: 
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Here Ty (My) is the y” total decay width given by 
PDG [29]. 

In previous analyses [37-40], the observed inclusive 
hadronic cross sections were fit with the theoretical one 
with the contributing terms in Eq. (8) to obtain the reso- 
nance parameters, which yield the 7)" cross section at the 
peak. The comparison of this cross section with the DD 
cross section measured by tagging D mesons yields the 
non-DD decay branching fraction. An inconsistency ex- 
ists in this procedure because in the fit to the total cross 
section, DD final states was assumed to saturate the y” 
decays. Since light hadrons have much lower thresholds 
than DD, a large non-DD branching fraction directly 
affects the shape of the resonance curve, and also indi- 
rectly through the energy-dependent width. Taking into 
account the non-DD decays, Eq. (9) should be revised 
by including another term, that is 


Ppo ae Pps LO = 
1+ (rppo) 1+ (rpps) pee | 3 


where Cyon-pp 1S proportional to the part of the width 
from non-DD decays, and 


Dyn (Myr) 
Cr = 3 3 
Ppo Ppt as C D| 
1+ (rpp)? 14+ (rpp+)? ROD i 
cm — p” 


With the Chon-pp term in the expression for Ty, the 
fitting of the resonance curve to extract the resonance 
parameters must be done together with the fitting of the 
DD cross section. In this procedure, the non-DD decay 
branching fraction is extracted together with the reso- 
nance parameters. 

However, this method subtracts DD cross section from 
the total inclusive one to get the non-DD cross section 
which is only a fraction of the total inclusive cross section. 
It suffers from the large uncertainties of the measure- 
ments, so it is difficult to obtain a statistically significant 
result. 


B. Exclusive method 


The calculations in section III show that those final 
states which are suppressed in ~' decays relative to J/w, 
and especially those ones enhanced in y’ decays, will 
show up in %” decays with maximum possible partial 
widths more than an order of magnitude greater than 
their widths in J/y or w’ decays. So these exclusive 
charmless modes should be searched for in 7” decays. 
This provides direct test of the calculations based on the 
2S-1D mixing scheme. Here we discuss three typical ex- 
clusive modes: VP mode which is suppressed in wy)’ de- 
cays relative to J/~, PP mode which is enhanced, and 
of particular interest, the BB and ¢f9(980) modes which 
observe the 12% rule. 

To measure the exclusive VP mode in w”" decays by 
ee~ experiments, the contribution from non-resonance 
virtual photon amplitude and its interference with the 
resonance must be treated with care. A recent study 
on the measurement of Y” + VP in ete experiments 
shows [24] that with the decay rate predicted by the S- 
and D-wave mixing, the interference between the three- 
gluon decay amplitude and the continuum one-photon 
amplitude leads to very small cross sections for some VP 
modes, e.g. pr and K**+ K- +<e.c., due to the destructive 
interference, but much larger cross sections for other VP 


+ 


modes, e.g. K “oK? +c.c. due to the constructive interfer- 
ence. In another word, although the branching fractions 
of po?r? and K “OR differ by only a fraction due to SU(3) 
symmetry breaking [41], their production cross sections 
in ete” collision differ by one to two orders of magnitude. 

Among the PP modes, there is the K2K® final state 
which decays only through strong interaction and does 
not couple to virtual photon [11, 41]. There is no compli- 
cation of electromagnetic interaction and the interference 
between it and the resonance. So the observed KLK? in 
ete experiment is completely from resonance decays. 
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In the 2-1D mixing scheme, with the BES recently re- 
ported KLK? branching fractions in J/2 [16] and y’ [17] 
decays as inputs, it is estimated [11] (1.2+0.7) x 10-® < 
Bi" + K&K®) < (3.8 + 1.1) x 1075 [42]. If there is no 
extra phase between (KK? |2°S1) and (K2K?|1°Dj), 
then its branching fraction is at the upper bound. With 
17.7 pb~! y” data, BES has set an upper limit [43], which 
is still beyond the sensitivity for testing the above pre- 
diction. More precise determination of this branching 
fraction is expected from the analysis based on larger 
data samples of CLEOc and BES-III. 

However, for other PP modes, or more generally other 
final states which are enhanced in 7' decays, in ete7 
experiments there is still the complication from the non- 
resonance virtual photon amplitude and its interference 
with the resonance. 

Of particular interest are the final states with Q(f) = 
1. These are the BB modes and the vector-scalar mode 
ofo(980) [44]. As discussed in subsection III C for Q(f) = 
1, there are two real and positive solutions with Rp = 
0.048 and Rr = 34.0. These two solutions are three 
orders of magnitude apart, their branching fractions are 
extremely sensitive to the relative phase between the 2S 
and 1D matrix elements. The BB branching fractions 
in J/w decays are at O(10~*), while ¢fo(980) is (3.2 + 
0.9) x 1074. If the physics solution of 7 is the larger one 
of the two real values, then the BB branching fraction 
in Y” decay would be at O(10~*) and ¢fo(980) would be 
4.2 x 107°, which can be observed in the 7” data sample 
over 1 fb7!. 


V. SUMMARY 


Based on the available experimental information of 
J/y and %' decays, we calculate the charmless decays 
of w" by virtue of the S- and D-wave charmonia mix- 
ing scheme which was proposed to explain the large 
yy" + ete partial width and the pr puzzle. We find 
that this leads to a possible large branching fraction, up 
to 13%, of the charmless final state in 7” decays. Al- 
though the calculation is semi-quantitative, it demon- 
strates that a large charmless branching fraction in y” 
decays can well be explained. 
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